Biochemical and transcriptomic evaluation of a 3D lung organoid platform for pre-clinical testing of active substances targeting senescence.

Chronic lung diseases such as chronic obstructive pulmonary disease and cystic fibrosis are incurable. Epithelial senescence, a state of dysfunctional cell cycle arrest, contributes to the progression of such diseases. Therefore, lung epithelial cells are a valuable target for therapeutic intervention. Here, we present a 3D airway lung organoid platform for the preclinical testing of active substances with regard to senescence, toxicity, and inflammation under standardized conditions in a 96 well format. Senescence was induced with doxorubicin and measured by activity of senescence associated galactosidase. Pharmaceutical compounds such as quercetin antagonized doxorubicin-induced senescence without compromising organoid integrity. Using single cell sequencing, we identified a subset of cells expressing senescence markers which was decreased by quercetin. Doxorubicin induced the expression of detoxification factors specifically in goblet cells independent of quercetin. In conclusion, our platform enables for the analysis of senescence-related processes and will allow the pre-selection of a wide range of compounds (e.g. natural products) in preclinical studies, thus reducing the need for animal testing.

Mutual Regulation of Transcriptomes between Murine Pneumocytes and Fibroblasts Mediates Alveolar Regeneration in Air-Liquid Interface Cultures.

Alveolar type 2 and club cells are part of the stem cell niche of the lung and their differentiation is required for pulmonary homeostasis and tissue regeneration. A disturbed crosstalk between fibroblasts and epithelial cells contributes to the loss of lung structure in chronic lung diseases. Therefore, it is important to understand how fibroblasts and lung epithelial cells interact during regeneration. Here, we analyzed the interaction of fibroblasts and the alveolar epithelium modeled in air-liquid interface cultures. Single-cell transcriptomics showed that cocultivation with fibroblasts leads to increased expression of type 2 markers in pneumocytes, activation of regulons associated with the maintenance of alveolar type 2 cells (e.g., Etv5), and transdifferentiation of club cells toward pneumocytes. This was accompanied by an intensified transepithelial barrier. Vice versa, the activation of NF-κB pathways and the CEBPB regulon and the expression of IL-6 and other differentiation factors (e.g., fibroblast growth factors) were increased in fibroblasts cocultured with epithelial cells. Recombinant IL-6 enhanced epithelial barrier formation. Therefore, in our coculture model, regulatory loops were identified by which lung epithelial cells mediate regeneration and differentiation of the alveolar epithelium in a cooperative manner with the mesenchymal compartment.

Ovulation is triggered by a cyclical modulation of gonadotropes into a hyperexcitable state.

Gonadotropes in the anterior pituitary gland are essential for fertility and provide a functional link between the brain and the gonads. To trigger ovulation, gonadotrope cells release massive amounts of luteinizing hormone (LH). The mechanism underlying this remains unclear. Here, we utilize a mouse model expressing a genetically encoded Ca2+ indicator exclusively in gonadotropes to dissect this mechanism in intact pituitaries. We demonstrate that female gonadotropes exclusively exhibit a state of hyperexcitability during the LH surge, resulting in spontaneous [Ca2+]i transients in these cells, which persist in the absence of any in vivo hormonal signals. L-type Ca2+ channels and transient receptor potential channel A1 (TRPA1) together with intracellular reactive oxygen species (ROS) levels ensure this state of hyperexcitability. Consistent with this, virus-assisted triple knockout of Trpa1 and L-type Ca2+ subunits in gonadotropes leads to vaginal closure in cycling females. Our data provide insight into molecular mechanisms required for ovulation and reproductive success in mammals.

Bitter taste cells in the ventricular walls of the murine brain regulate glucose homeostasis.

The median eminence (ME) is a circumventricular organ at the base of the brain that controls body homeostasis. Tanycytes are its specialized glial cells that constitute the ventricular walls and regulate different physiological states, however individual signaling pathways in these cells are incompletely understood. Here, we identify a functional tanycyte subpopulation that expresses key taste transduction genes including bitter taste receptors, the G protein gustducin and the gustatory ion channel TRPM5 (M5). M5 tanycytes have access to blood-borne cues via processes extended towards diaphragmed endothelial fenestrations in the ME and mediate bidirectional communication between the cerebrospinal fluid and blood. This subpopulation responds to metabolic signals including leptin and other hormonal cues and is transcriptionally reprogrammed upon fasting. Acute M5 tanycyte activation induces insulin secretion and acute diphtheria toxin-mediated M5 tanycyte depletion results in impaired glucose tolerance in diet-induced obese mice. We provide a cellular and molecular framework that defines how bitter taste cells in the ME integrate chemosensation with metabolism.

Intra-pituitary follicle-stimulating hormone signaling regulates hepatic lipid metabolism in mice.

Inter-organ communication is a major hallmark of health and is often orchestrated by hormones released by the anterior pituitary gland. Pituitary gonadotropes secrete follicle-stimulating hormone (FSH) and luteinizing hormone (LH) to regulate gonadal function and control fertility. Whether FSH and LH also act on organs other than the gonads is debated. Here, we find that gonadotrope depletion in adult female mice triggers profound hypogonadism, obesity, glucose intolerance, fatty liver, and bone loss. The absence of sex steroids precipitates these phenotypes, with the notable exception of fatty liver, which results from ovary-independent actions of FSH. We uncover paracrine FSH action on pituitary corticotropes as a mechanism to restrain the production of corticosterone and prevent hepatic steatosis. Our data demonstrate that functional communication of two distinct hormone-secreting cell populations in the pituitary regulates hepatic lipid metabolism.

Transcriptional and Epigenetic Consequences of DMSO Treatment on HepaRG Cells.

Dimethyl sulfoxide (DMSO) is used to sustain or favor hepatocyte differentiation in vitro. Thus, DMSO is used in the differentiation protocol of the HepaRG cells that present the closest drug-metabolizing enzyme activities to primary human hepatocytes in culture. The aim of our study is to clarify its influence on liver-specific gene expression. For that purpose, we performed a large-scale analysis (gene expression and histone modification) to determine the global role of DMSO exposure during the differentiation process of the HepaRG cells. The addition of DMSO drives the upregulation of genes mainly regulated by PXR and PPARα whereas genes not affected by this addition are regulated by HNF1α, HNF4α, and PPARα. DMSO-differentiated-HepaRG cells show a differential expression for genes regulated by histone acetylation, while differentiated-HepaRG cells without DMSO show gene signatures associated with histone deacetylases. In addition, we observed an interplay between cytoskeleton organization and EMC remodeling with hepatocyte maturation.

Identification of an FXR-modulated liver-intestine hybrid state in iPSC-derived hepatocyte-like cells.

BACKGROUND & AIMS: Pluripotent stem cell (PSC)-derived hepatocyte-like cells (HLC) have enormous potential as a replacement for primary hepatocytes in drug screening, toxicology and cell replacement therapy, but their genome-wide expression patterns differ strongly from primary human hepatocytes (PHH). METHODS: We differentiated human induced pluripotent stem cells (hiPSC) via definitive endoderm to HLC and characterized the cells by single-cell and bulk RNA-seq, with complementary epigenetic analyses. We then compared HLC to PHH and publicly available data on human fetal hepatocytes (FH) ex vivo; we performed bioinformatics-guided interventions to improve HLC differentiation via lentiviral transduction of the nuclear receptor FXR and agonist exposure. RESULTS: Single-cell RNA-seq revealed that transcriptomes of individual HLC display a hybrid state, where hepatocyte-associated genes are expressed in concert with genes that are not expressed in PHH - mostly intestinal genes - within the same cell. Bulk-level overrepresentation analysis, as well as regulon analysis at the single-cell level, identified sets of regulatory factors discriminating HLC, FH, and PHH, hinting at a central role for the nuclear receptor FXR in the functional maturation of HLC. Combined FXR expression plus agonist exposure enhanced the expression of hepatocyte-associated genes and increased the ability of bile canalicular secretion as well as lipid droplet formation, thereby increasing HLCs' similarity to PHH. The undesired non-liver gene expression was reproducibly decreased, although only by a moderate degree. CONCLUSION: In contrast to physiological hepatocyte precursor cells and mature hepatocytes, HLC co-express liver and hybrid genes in the same cell. Targeted modification of the FXR gene regulatory network improves their differentiation by suppressing intestinal traits whilst inducing hepatocyte features. LAY SUMMARY: Generation of human hepatocytes from stem cells represents an active research field but its success is hampered by the fact that the stem cell-derived 'hepatocytes' still show major differences to hepatocytes obtained from a liver. Here, we identified an important reason for the difference, specifically that the stem cell-derived 'hepatocyte' represents a hybrid cell with features of hepatocytes and intestinal cells. We show that a specific protein (FXR) suppresses intestinal and induces liver features, thus bringing the stem cell-derived cells closer to hepatocytes derived from human livers.

Transmission of trained immunity and heterologous resistance to infections across generations.

Intergenerational inheritance of immune traits linked to epigenetic modifications has been demonstrated in plants and invertebrates. Here we provide evidence for transmission of trained immunity across generations to murine progeny that survived a sublethal systemic infection with Candida albicans or a zymosan challenge. The progeny of trained mice exhibited cellular, developmental, transcriptional and epigenetic changes associated with the bone marrow-resident myeloid effector and progenitor cell compartment. Moreover, the progeny of trained mice showed enhanced responsiveness to endotoxin challenge, alongside improved protection against systemic heterologous Escherichia coli and Listeria monocytogenes infections. Sperm DNA of parental male mice intravenously infected with the fungus C. albicans showed DNA methylation differences linked to immune gene loci. These results provide evidence for inheritance of trained immunity in mammals, enhancing protection against infections.

Prolactin-sensitive olfactory sensory neurons regulate male preference in female mice by modulating responses to chemosensory cues.

Chemosensory cues detected in the nose need to be integrated with the hormonal status to trigger appropriate behaviors, but the neural circuits linking the olfactory and the endocrine system are insufficiently understood. Here, we characterize olfactory sensory neurons in the murine nose that respond to the pituitary hormone prolactin. Deletion of prolactin receptor in these cells results in impaired detection of social odors and blunts male preference in females. The prolactin-responsive olfactory sensory neurons exhibit a distinctive projection pattern to the brain that is similar across different individuals and express a limited subset of chemosensory receptors. Prolactin modulates the responses within these neurons to discrete chemosensory cues contained in male urine, providing a mechanism by which the hormonal status can be directly linked with distinct olfactory cues to generate appropriate behavioral responses.

Rapid Base-Specific Calling of SARS-CoV-2 Variants of Concern Using Combined RT-PCR Melting Curve Screening and SIRPH Technology.

BACKGROUND: The emergence of novel variants of concern of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demands fast and reliable detection of such variants in local populations. METHODS: Here we present a cost-efficient and fast workflow combining a prescreening of SARS-CoV-2-positive samples using reverse transcription polymerase chain reaction melting curve analysis with multiplexed IP-RP-HPLC-based single nucleotide primer extensions. RESULTS: The entire workflow from positive SARS-CoV-2 testing to base-specific identification of variants requires about 24 hours. CONCLUSIONS: We applied the sensitive method to monitor local variant of concern outbreaks in SARS-CoV-2-positive samples collected in a confined region of Germany.

CpG content-dependent associations between transcription factors and histone modifications.

Understanding the factors that underlie the epigenetic regulation of genes is crucial to understand the gene regulatory machinery as a whole. Several experimental and computational studies examined the relationship between different factors involved. Here we investigate the relationship between transcription factors (TFs) and histone modifications (HMs), based on ChIP-seq data in cell lines. As it was shown that gene regulation by TFs differs depending on the CpG class of a promoter, we study the impact of the CpG content in promoters on the associations between TFs and HMs. We suggest an approach based on sparse linear regression models to infer associations between TFs and HMs with respect to CpG content. A study of the partial correlation of HMs for the two classes of high and low CpG content reveals possible CpG dependence and potential candidates for confounding factors in our models. We show that the models are accurate, inferred associations reflect known biological relationships, and we give new insight into associations with respect to CpG content. Moreover, analysis of a ChIP-seq dataset in HepG2 cells of the HM H3K122ac, an HM about little is known, reveals novel TF associations and supports a previously established link to active transcription.

Prediction of single-cell gene expression for transcription factor analysis.

BACKGROUND: Single-cell RNA sequencing is a powerful technology to discover new cell types and study biological processes in complex biological samples. A current challenge is to predict transcription factor (TF) regulation from single-cell RNA data. RESULTS: Here, we propose a novel approach for predicting gene expression at the single-cell level using cis-regulatory motifs, as well as epigenetic features. We designed a tree-guided multi-task learning framework that considers each cell as a task. Through this framework we were able to explain the single-cell gene expression values using either TF binding affinities or TF ChIP-seq data measured at specific genomic regions. TFs identified using these models could be validated by the literature. CONCLUSION: Our proposed method allows us to identify distinct TFs that show cell type-specific regulation. This approach is not limited to TFs but can use any type of data that can potentially be used in explaining gene expression at the single-cell level to study factors that drive differentiation or show abnormal regulation in disease. The implementation of our workflow can be accessed under an MIT license via https://github.com/SchulzLab/Triangulate.

Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data.

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-seq, ATAC-seq and NOMe-seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However, we also find a high proportion of assay specific NDRs that are often 'called' by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-seq and DNase I-seq provide a superb high NDR calling rate for relatively low sequencing costs in comparison to NOMe-seq, NOMe-seq singles out for its genome-wide coverage allowing to not only detect NDRs but also endogenous DNA methylation and as we show here genome wide segmentation into heterochromatic B domains and local phasing of nucleosomes outside of NDRs. In summary, our comparisons strongly suggest to consider assay specific differences for the experimental design and for generalized and comparative functional interpretations.

GABA is a modulator, rather than a classical transmitter, in the medial nucleus of the trapezoid body-lateral superior olive sound localization circuit.

KEY POINTS: The lateral superior olive (LSO), a brainstem hub involved in sound localization, integrates excitatory and inhibitory inputs from the ipsilateral and the contralateral ear, respectively. In gerbils and rats, inhibition to the LSO reportedly shifts from GABAergic to glycinergic within the first three postnatal weeks. Surprisingly, we found no evidence for synaptic GABA signalling during this time window in mouse LSO principal neurons. However, we found that presynaptic GABAB Rs modulate Ca2+ influx into medial nucleus of the trapezoid body axon terminals, resulting in reduced synaptic strength. Moreover, GABA elicited strong responses in LSO neurons that were mediated by extrasynaptic GABAA Rs. RNA sequencing revealed highly abundant δ subunits, which are characteristic of extrasynaptic receptors. Whereas GABA increased the excitability of neonatal LSO neurons, it reduced the excitability around hearing onset. Collectively, GABA appears to control the excitability of mouse LSO neurons via extrasynaptic and presynaptic signalling. Thus, GABA acts as a modulator, rather than as a classical transmitter. ABSTRACT: GABA and glycine mediate fast inhibitory neurotransmission and are coreleased at several synapse types. Here we assessed the contribution of GABA and glycine in synaptic transmission between the medial nucleus of the trapezoid body (MNTB) and the lateral superior olive (LSO), two nuclei involved in sound localization. Whole-cell patch-clamp experiments in acute mouse brainstem slices at postnatal days (P) 4 and 11 during pharmacological blockade of GABAA receptors (GABAA Rs) and/or glycine receptors demonstrated no GABAergic synaptic component on LSO principal neurons. A GABAergic component was absent in evoked inhibitory postsynaptic currents and miniature events. Coimmunofluorescence experiments revealed no codistribution of the presynaptic GABAergic marker GAD65/67 with gephyrin, a postsynaptic marker for GABAA Rs, corroborating the conclusion that GABA does not act synaptically in the mouse LSO. Imaging experiments revealed reduced Ca2+ influx into MNTB axon terminals following activation of presynaptic GABAB Rs. GABAB R activation reduced the synaptic strength at P4 and P11. GABA appears to act on extrasynaptic GABAA Rs as demonstrated by application of 4,5,6,7-tetrahydroisoxazolo[5,4-c]pyridin-3-ol, a δ-subunit-specific GABAA R agonist. RNA sequencing showed high mRNA levels for the δ-subunit in the LSO. Moreover, GABA transporters GAT-1 and GAT-3 appear to control extracellular GABA. Finally, we show an age-dependent effect of GABA on the excitability of LSO neurons. Whereas tonic GABA increased the excitability at P4, leading to spike facilitation, it decreased the excitability at P11 via shunting inhibition through extrasynaptic GABAA Rs. Taken together, we demonstrate a modulatory role of GABA in the murine LSO, rather than a function as a classical synaptic transmitter.

The role of TCF3 as potential master regulator in blastemal Wilms tumors.

Wilms tumors are the most common type of pediatric kidney tumors. While the overall prognosis for patients is favorable, especially tumors that exhibit a blastemal subtype after preoperative chemotherapy have a poor prognosis. For an improved risk assessment and therapy stratification, it is essential to identify the driving factors that are distinctive for this aggressive subtype. In our study, we compared gene expression profiles of 33 tumor biopsies (17 blastemal and 16 other tumors) after neoadjuvant chemotherapy. The analysis of this dataset using the Regulator Gene Association Enrichment algorithm successfully identified several biomarkers and associated molecular mechanisms that distinguish between blastemal and nonblastemal Wilms tumors. Specifically, regulators involved in embryonic development and epigenetic processes like chromatin remodeling and histone modification play an essential role in blastemal tumors. In this context, we especially identified TCF3 as the central regulatory element. Furthermore, the comparison of ChIP-Seq data of Wilms tumor cell cultures from a blastemal mouse xenograft and a stromal tumor provided further evidence that the chromatin states of blastemal cells share characteristics with embryonic stem cells that are not present in the stromal tumor cell line. These stem-cell like characteristics could potentially add to the increased malignancy and chemoresistance of the blastemal subtype. Along with TCF3, we detected several additional biomarkers that are distinctive for blastemal Wilms tumors after neoadjuvant chemotherapy and that may provide leads for new therapeutic regimens.

Integrative analysis of single-cell expression data reveals distinct regulatory states in bidirectional promoters.

BACKGROUND: Bidirectional promoters (BPs) are prevalent in eukaryotic genomes. However, it is poorly understood how the cell integrates different epigenomic information, such as transcription factor (TF) binding and chromatin marks, to drive gene expression at BPs. Single-cell sequencing technologies are revolutionizing the field of genome biology. Therefore, this study focuses on the integration of single-cell RNA-seq data with bulk ChIP-seq and other epigenetics data, for which single-cell technologies are not yet established, in the context of BPs. RESULTS: We performed integrative analyses of novel human single-cell RNA-seq (scRNA-seq) data with bulk ChIP-seq and other epigenetics data. scRNA-seq data revealed distinct transcription states of BPs that were previously not recognized. We find associations between these transcription states to distinct patterns in structural gene features, DNA accessibility, histone modification, DNA methylation and TF binding profiles. CONCLUSIONS: Our results suggest that a complex interplay of all of these elements is required to achieve BP-specific transcriptional output in this specialized promoter configuration. Further, our study implies that novel statistical methods can be developed to deconvolute masked subpopulations of cells measured with different bulk epigenomic assays using scRNA-seq data.

Epigenomic map of human liver reveals principles of zonated morphogenic and metabolic control.

A deeper epigenomic understanding of spatial organization of cells in human tissues is an important challenge. Here we report the first combined positional analysis of transcriptomes and methylomes across three micro-dissected zones (pericentral, intermediate and periportal) of human liver. We identify pronounced anti-correlated transcriptional and methylation gradients including a core of 271 genes controlling zonated metabolic and morphogen networks and observe a prominent porto-central gradient of DNA methylation at binding sites of 46 transcription factors. The gradient includes an epigenetic and transcriptional Wnt signature supporting the concept of a pericentral hepatocyte regeneration pathway under steady-state conditions. While donors with non-alcoholic fatty liver disease show consistent gene expression differences corresponding to the severity of the disease across all zones, the relative zonated gene expression and DNA methylation patterns remain unchanged. Overall our data provide a wealth of new positional insights into zonal networks controlled by epigenetic and transcriptional gradients in human liver.

A comprehensive analysis of 195 DNA methylomes reveals shared and cell-specific features of partially methylated domains.

BACKGROUND: Partially methylated domains are extended regions in the genome exhibiting a reduced average DNA methylation level. They cover gene-poor and transcriptionally inactive regions and tend to be heterochromatic. We present a comprehensive comparative analysis of partially methylated domains in human and mouse cells, to identify structural and functional features associated with them. RESULTS: Partially methylated domains are present in up to 75% of the genome in human and mouse cells irrespective of their tissue or cell origin. Each cell type has a distinct set of partially methylated domains, and genes expressed in such domains show a strong cell type effect. The methylation level varies between cell types with a more pronounced effect in differentiating and replicating cells. The lowest level of methylation is observed in highly proliferating and immortal cancer cell lines. A decrease of DNA methylation within partially methylated domains tends to be linked to an increase in heterochromatic histone marks and a decrease of gene expression. Characteristic combinations of heterochromatic signatures in partially methylated domains are linked to domains of early and middle S-phase and late S-G2 phases of DNA replication. CONCLUSIONS: Partially methylated domains are prominent signatures of long-range epigenomic organization. Integrative analysis identifies them as important general, lineage- and cell type-specific topological features. Changes in partially methylated domains are hallmarks of cell differentiation, with decreased methylation levels and increased heterochromatic marks being linked to enhanced cell proliferation. In combination with broad histone marks, partially methylated domains demarcate distinct domains of late DNA replication.